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ABSTRACT 



The ratio cf geostrophic drag coefficients C/C a under 
various stability conditions is obtained as a function of 
Richardson number. A drag coefficient Cf is defined in 
relation to the mechanical mixing length, and its ratio to 
the actual drag coefficient C is related to the Richardson 
number. By using observational recent data on the normal- 
ized logarithmic wind shear of Monin-Obukhov one can obtain 
Cf/C a . 

The writer is deeply indebted to Dr. P. L» Martin 
(Professor of Meteorology) for his suggestion of the topic 
and his continued help throughout the investigation and 
during the preparation of this paper. 
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1. Introduction. 



One of the important problems in meteorology is to 
establish relationships between the turbulent characteris- 
tics of the lower atmosphere and the large-scale synoptic 
parameters. A significant turbulent characteristic is the 
ground drag, which enables one to specify the wind profile, 
the eddy diffusivity and the local energy dissipation in 
the surface layer, by knowing the topographical characteris- 
tics i.e., the roughness parameter of a particular place. 

Lettau [ 7 ] had made an attempt to relate the drag 
coefficients with the static stability represented by a 
dimensionless parameter called Richardson number. His 
investi ations were based on 1953 O'Neill, Nebraska, data. 

It is the main purpose of this paper to verify and extend 
his investigations using 1956 O'Neill data. It is important 
to note that this wo rk deals with a relatively smooth 
surface, therefore the "zero-plane displacement" which is 
sometimes used in the adiabatic-wind profile has been 
neglected, a procedure which greatly simplifies the com- 
putations. 
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2. The Nature of the Data. 



The data used in this study were obtained by the 
Project Prairie Grass ^l] conducted near O'Neill, Nebraska 
durinp the summer 1956. Since this study rakes use of some 
of the published data, i.e., wind speed, temperature and 
surface maps, it is necessary to describe some aspects of 
the instrumentation. The instrumentation employed was that 
of the mobile rni crorneteorological station of the Texas A&M 
group. The station has a slender aluminum mast supporting 
six anemometers at heights of 8, 4, 2, 1, 0.5 and 0.25 meter 
A similar mast supported seven temperature-measuring, radia- 
tion-shielded copper-constantan thermocouple junctions at 
heights 8, 4, 2, 1, 0.5, 0.25 and 0.125 meters. In addition 
the published report of Project Prairie Grass includes 
sectional sea level pressure maps at times nearly synoptic 
with the micrometeorological data. These maps revealed some 
small-scale features of the circulation near O'Neill which 
did not appear on the larger-scale facsimile maps. The micro 
meteorological data periods were selected to be nearly simul 
taneous with the available maps. At times, linear interpola 
tion between the available charts was used for computations 
of geostrophic wind speeds. 

Regarding the errors of measurement, two' types of tem- 
perature error occur: i.e., calibration and radiation error. 
An estimate for calibration error was 0.05OC. The estimated 
radiation error was 0.1°C with the understanding tl at, in 
daytime with clear skies, and low windspeed, all measured 



2 



air temperatures overestimate the true temperature. At night 
with clear skies the measured temperatures would he lower 
than the ambient temperature. For windspeed, measurements 
were dependent upon different techniques of calibrations 
according to the windforce. All measurements used in this 
paper were twenty-minute means. The reader is referred to 
[l] for more details concerning the instrumentation. 



5. Background Theory. 

Extensive use will be made of Prandtl's turbulence 
model, because of its relative simplicity. This model is 
described in many standard texts. According to this theory, 
the f ollowin major results emerge: 



a = X' # 

dz 



x' = - p 



(1) 



( 2 ) 



where JUL' , fuh‘ , are the turbulent velocity fluctuations in 
the direction of the mean wind and in the vertical, respec- 
tively. £ is the eddy mixing" length, while & is the root- 
mean-souare mixing len th. The eddy stress X' is then 
considered to he very nearly a constant X c within a layer 
called the surface laj^er. 

In the case of the fully-rough, neutral surface layer, 
Rossby has shown thatXski, so that (2) may he integrated 
to give 

U.* = /tvu 5. (3) 



where U,* & •= VT Jf is the friction velocity in a neutral 
surface layer, Z 0 is the rou-hness parameter and k=0.38 



is the Von Karman constant. 

In the non-neutral surface layer, numerous theoretical 
profiles have been derived. However Deacon has proposed the 
semi-empirical relationship 



da. __ ul» ( iX* 

where the friction velocity u.*. is constant in the 



(4) 



non-neutral surface layer. The parameter 6 is a decreasing 
function of Richardson number. Actually Davidson and 
larad [ 4 ] have shown that |0-l\ increases somewhat with 
heiht. However (4) still affords a useful tool as an 
overall or bulk-relationship. Two useful res'. Its may be 
derived from (4). These are: 

(i) the integrated wind profile; 

i-ft “i 



U - 



U., 



[(*)- 



(5) 

kO-G) 

(ii) the value JU.* by the limitin process indicated 
in eouation (7) below 



a. 



I, X- 6 u= / 6 (U 

- * Z-*Zo“ WW* 



(5) 



so that 



V |ji- 

* 2-*4 0(W) 



(7) 



Actually the lowest level at which wind data were available 
was Z s 25 cm and Z D is the rouphness ja rameter determined 
under the time-nearest neutral conditions. U.* may therefore 
be approximated by the finite-difference expression 
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XX* — 



0. 38 XXjy 



( 6 ) 



JUa5- jUz { 

Following a procedure similar to that employed by David 
son and Farad [4] in neutral cases, it may be shown that 

jL- ”' Li; *(*»*) 

where 



- 






(9) 



I el 



A(bvz) = u. + U + U ago 



( 10 ) 



and 



ALL = (u 6 -u.,) + (a s .-u 2 ) + ( u . ii .a^ (id 

for mean-wind measurements at 25, 50, 100, 200, 400, 800 cms 



4. Oeostrophic Drag Coefficients 0 

Lettau M v as defined the peostrophic drag coefficient 
as » 

C = a, /v 3 (is) 

For neutral stability, the corresponding drag coefficient 
is defined as 

0^= XA.*o./V 3 (13) 

Neutral wind profiles at O'Neill may be obtained by finding 
those cases for which the potential temperature profiles 
we characterized byds/dzuO in the layer 25 cm - 2C0 cm. 
Altogether 15 such profiles were found. Assuming a lo'arith 
mic wind profile can be fitted to the wind data, * 1 * 0 . of 

equation (3) may be obtained as 

icAU. L ‘ “ 



U *' J - = A(jUvz) 



VhU. 



(14) 



By analogy with flow in circular conduits, Lettau [y} 
has suggested a drag-coefficient Ca for neutral cases of 
th e f o rm 






Q.IQ4 



(15) 



“ jLo 3 (C^R 0 )«2.24 

where Ro may he called the "surface Rossby number" defined 

as the non-dimensional expression 



Ko=V 3 /f Zo , V 3 ,*±£ (1«) 

The constants of eouation (15) were obtained by Lettau by 

least-square methods applied to a wide variety of neutral- 

wind profiles for different sites. 

Lettau £ 7 } goes on to give values of C/C a for different 

values R1 ]_qo> the Richardson number at 100 cm, for the 1953 

O'Neill data. He obtains a curve as shown in fig. 1. 




.Ratio of geostrophic drag coefficient for non-adiabati c and 
adiabatic conditions as a function of Richardson number at 
Z= 100 cm, at O'Neill, Nebraska 1953 (After Lettau £7} ). 

It must be pointed out that C Q here represents the geostro- 
phic drag coefficient for any surface layer of drag coeffi- 
cient C, if this layer were suddenly to be converted into 



a neutral layer without changing Vg/fz 0 . 

A primary objective of this paper was to test eauation 
(15) as a means of determining C a based on the 1956 O'Neill 
data. Also application of non-neutral orofile data to eoua- 

v 

tions (7) and (12) afforded values of C under identical 
roughness conditions. The comparison of Lettau's 1953 graph 
of C/C a versus Riioo with that obtained by the writer at 
Hll4 0 (to be more exact, the Richardson number at 141 0 4cm) 
based on the 1956 O'Neill data is shown in fig. 3. 

If one knows a value of Ro = Vg/fz 0 , it is possible to 
solve eauation (15) for C a by numerical methods. It is 
shown in the Appendix that a numerical solution C a of eaua- 
tion (15) may be obtained iteratively from the eauation 

y.* (a. 1026 iUy ; + y- + - 2.3o3$-2 c <] fte ) , 

v - y, - — — : ‘ — — — - (17) 

• 1-4- y* 

where i-th iterant to 0.2395/C a . The first estimate 
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The chief inaccuracy in this approach appears to be uncer- 
tainty in the numerical constants appearing in eauation(l5) , 
especially at low geostrophic wind speeds. Another source 
of error is the necessity in some cases, of interpolating 
Vg between two maps 6 or 12 hours apart in order to get the 
estimated value of Vg at the time of micrometeorolopical 
data-run. 

A second approach to determine the relationship between 
C and C a was used. Based on 15 neutral cases at O'Neill, 
Nebraska, values of C a were computed using equations (13) 
and (14). Values of C a were then plotted versus log Ro, 
and the resulting best-fitting curve compared to Lettau's, 
which is replotted on the same graph. The comparison will 
be discussed in later section. 

The Richardson number has been used very extensively 
as a measure of stability in the micro-meteorology. Its 
definition is 






j 2>e/d£ 



(18) 



6 ( du./d2) a 

This number is proportional to the ratio of work done against 
static stability to the work done by eddy stress. 

For instruments arranged at eoual logarithmic spacin? , 
as in the 1953 and 1956 O'Neill data, a convenient finite- 
difference form of Ri at Z = has been proposed by 

Lettau [sj 

r-\ / A a ? . . ( R. - Si JLl, 9. 

(19) 



Kl1 J= n&a-V 



This value is understood to be applicable at the geometric 
of Z]_ and Z 2 , namely Z\, 2 ~ *^"ziZ2 . ©2 and 0-2 in the 
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mean 



1956 data are the 20-minute mean potential temperature and 
windspeed at level 2, and an analogous meaning applies for 
symbols with subscript one. 

Finally, Monin and Obukhov M have introduced a use- 
ful non-dimensional parameter S defined by 

_ kz (g 0) 

iz. 

Note that the numerator of S has the form 

k*dCL/dz - (21) 

which, according to the Prandtl model, may be regarded as 
representing as here defined represents that 

part of tX* due to turbulent fluctuations associated with 
the mechanical turbulent mixing lenpth 1^-ski. The expression 
for 3 may then be written 

S 

The drap coefficient defined by Cf -z. U.*^/ Vg is not one 
that can be measured directly, but is one whose value can be 
inferred by empirical information regarding S. This is 
dealt with in more detail in section 8. This interpretation 
of as kz appears to have been piven first by 

Bu singer [sj . 

5. Comparison of C/C a for 1952 and 1956 data. 

In computin'" C a , Lettau’3 empirical formula, ecuation 
(17) was first used with the 1956 data. It was necessary 
to compute Z 0 , usinp ecuation (9) from the time-nearest 
neutral wind profiles and Vp from interpolated surface 



_ / Vg _ 

~ u»/v / 3 - c 



( 22 ) 




9 



cherts. In computing U.*and C, eouation (8) and (12) were 
used. The value of' ftip 4 Q is also included using equation 
(19). The results for all data-days employed are given in 
table 1, and the groupings ol Rii40 in classes are given in 
table 2. The graph of C/Ca versus Rii40 Is given in fig. 3, 
with Lettau’s curve superimposed. The agreement with 

Lettau’s results is not as good as might have been expected. 
For instance the value C/C a - 1.75 was found at Ri]_4Q-0 
contrasted with Lettau's value C/C a -1.00 at RipQO”® 0 More- 
over the maximum value of C/C a found in this paper exceeds 
Lettau’s maximum by a comparable factor but with Rii40 - 
-0.055, compared with RipoO- ~0»007, according to Lettau. 

Possible sources of error lie in the values of Vg, 
which affect the values of C and C a given by eouations (12) 
and (15) respectively, as well as in the value of Z 0 . The 
values of Z 0 determined here ranged from 0.35cm to 3.1 cm, 
as contrasted with values quoted by Blackadar et al of 

Z 0 — 0.6 cm for profiles in approximately the same data- 
period. 
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TABLE I - Values of C/C a 
Date CST Z Q (cms) U*fcrn/sec) 
10 Jul 1955 



1505 




1.10 


24.67 


1405 




1.10 


35.85 


1505 




1.10 


31.24 


1530 




1.10 


30.99 


1605 




1.10 


35.75 


1705 




1.10 


35.00 


1905 




1.10 


44.00 


2105 




1.10 


21.03 


Jul 


1956 






0930 




3.14 


49.26 


1105 




3.14 


42.85 


1205 




3.14 


46.69 


1305 




3.14 


52.92 


1605 




3.14 


52.55 


1705 




3.14 


59 . 69 


1805 




3.14 


67.38 


2130 




3.14 


31.03 


2305 




3.14 


31.86 


i Jul 


1956 






0205 




1.49 


2.15 


0405 




1.49 


12.66 


0805 




1.49 


50.49 


0905 




1.49 


52.24 


1230 




0.70 


58.87 


1405 




0.70 


65 24 


1605 




0.70 


62.80 


2330 




0.70 


48.72 


Aup 


1S56 






1805 




0.35 


29.53 


2105 




0.35 


11.65 


2205 




0.35 


25.00 



eauations ( 3 ), ( 12 ), ( 15 ). 



( kts ) 


C 


C a 


C / C a ^140 


60 


0.036 


0.034 


1.08 - 0.46 


,80 


0 . 055 


0.034 


1.65 - 0.17 


00 


0.052 


0.034 


1.52 - 0.28 


64 


0.040 


0.033 


1.20 - 0.25 


60 


0.053 


0.034 


1.54 - 0.20 


80 


0.065 


0.035 


1.87 - 0.05 


52 


0.094 


0.035 


2.41 - 0 . C2 


00 


0.035 


0.034 


1.02 C. 06 



44 


0.051 


0.037 


1.38 


- 0.15 


00 


0.045 


0.037 


1.22 


- 0.21 


00 


0.052 


0.037 


1.39 


- 0.16 


00 


0.062 


0.038 


1.66 


- 0.09 


20 


0.065 


0.038 


1.72 


- 0.14 


00 


0.075 


0.038 


1.97 


- 0.06 


80 


0.085 


0.038 


2.25 


- 0.02 


52 


0.059 


0 . C40 


1.47 


0.04 


52 


0.085 


0.042 


2. Cl 


0.03 



80 


0.009 


0.040 


0.22 


0.10 


60 


0.026 


0.037 


0.72 


0.05 


00 


0.042 


0.034 


1.22 


- 0.10 


00 


0.040 


0.032 


1.24 


- 0.01 


00 


0.038 


0.029 


1.27 


- 0.02 


00 


0.028 


0.028 


1.02 


- 0.02 


00 


0.023 


0.027 


0.85 


- 0.01 


80 


0.022 


0.028 


0.78 


0.00 



60 


0.056 


0.030 


1.84 


0.00 


80 


0.022 


0.030 


0.71 


0.07 


00 


0.045 


0.030 


1 . 50 


0.01 



by 

Vf 

15 

12 

12 

15 

11 

10 

10 

12 

19 

19 

18 

17 

16 

16 

15 

9 

7 

4 

9 

24 

26 

52 

46 

54 

44 

10 

10 

11 
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TABL8 1 - (Cont'd) 



Date CST Z 0 (cms) ujpm/sec) 
7 Aup 1956 





0405 


3.06 


57.72 




0805 


3.06 


63.14 




0905 


3.06 


109.10 




1105 


3.06 


69.66 




1530 


0, 62 


59.79 




1605 


0.62 


64.83 




1805 


0.62 


36.02 




2205 


0. 62 


23.87 




2205 


0.62 


9.47 


27 


Auf 1956 








1230 


0.55 


42.96 




1505 


0.55 


38.09 




1705 


0.55 


40.08 




1930 


0.55 


1.79 




2205 


0.55 


2.98 


29 


Aug 1956 








2135 


1.53 


24.09 



Vg (kts) 


C 


c a 


C/C a ill 140 


11.88 


0.097 


0.039 


2.47 


0.01 


12.88 


0.098 


0.039 


2.52 -0.04 


13.08 


0.167 


0.039 


4.30 -0.02 


13.48 


0.103 


0.039 


2.67 -0.08 


15.84 


0.075 


0.031 


2.44 -0.04 


15.00 


0.086 


0.031 


2.77 -0.01 


15.76 


0.046 


0.031 


1.47 =0.01 


14.00 


0.034 


0.031 


1.09 


0.05 


12.40 


0.015 


0.032 


0.48 


0.13 



26.48 


0.032 


0.029 


1.12 


-0.07 


24.80 


0.031 


0.029 


1.06 


-0.02 


25.20 


0.032 


0.029 


1.10 


-0.02 


15.92 


0.002 


0.031 


0.07 


0.15 


26. 60 


0.002 


0.029 


0.07 


0.04 


7.20 


0.057 


CO 

to 

o 

0 

o 


1.75 


0.03 



TABLS 2 - Class means of Ri ]_40 vs C/C a 



Ri 140 


- 0.43 


- 0.37 


- 0.31 


- 0,25 


- 0.19 


- 0,12 


- 0.06 


0.00 


0.06 


0.13 


c/C a 


1.08 


- 


1.52 


1.20 


1.45 


1.44 


2.04 


1.72 


1.11 


0.26 
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Pig. 3 

Ratio of geostrophic drag coefficient for non-adiabatic and 
adiabatic conditions as a function of Richardson number 
at Z = 140 cm, at O’Neill, Nebraska 1956. 
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6. Revised formula for C a 

Hitherto C a has been computed by using Lettau’s empirical 
function for C a , eouation (15). Since Lettau himself realized 
that the numerical constants in (15) were not very exact 7 : , 
an attempt is made here to obtain an improved relationship 
between C a and Ro by usinr only 1956 O’Neill data. The first 
step again consists in computing Z 0 for the time-nearest 
reutral wind profile according to eouation (9) 0 As before,, 
U.**. results from equation (14) with ks 0 o 38, so that 



— 0.06 J AUlc 

Eouation (13) is then applied to determine C a „ The results 
of the computations are shown in table 3. 

An empirical function for C a was then obtained by the 
followin'" procedure: 

Let log Ro-x and C a ^ y and suppose the desired function has 
the form 

y -be (24) 

where y means the desired function, while b and <g are cop- 
stants to be determined. Taking the log of (24) one obtains 



mita-b+cx (25) 

To simplify the notation, let SLy.y and IU sa, and there- 
fore the desired function becomes 

Y = Oi + C)!. (26) 

Using the least-squares method, the normal ecuations are 

Y — CL— C-X- - O 

2Y«-aZ*, -cK.o (2,) 

... is» 

Substituting the values of observations (table 3) in the 
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normal equations one obtains 



3.51 - a - 7 „ 21 c = 0 



38190-103. 28a - 785.18c =0 



from which, it follows that 



a = -0.85 or b - 0.427 



c = -0.369 



Substitutinf these constants in (24), and further by chanp~ 
inf y to C a , x to lop Ho, the expression for C a can be 
obtained 



The fraph of eouation (28) is shown in fif. 4. compared with 
Lettau’s fraph. From fif. 4 it is evident that the revised 
formula for C a tends to underestimate the draf at larger 
values of liossby number, and overestimate it at smaller xio 
(lof ito<7.1). It should be noted, however, that there was 
considerable scatter of the observations relative to the line 
of best fit. 

If one assumes that eouation (28) is a universal formula 
for eostrophic drap coefficient in riddle and higher latitudes 
under neutral conditions, one can then write 



with C a fiven by (28). Hence, adiabatic wind speed at level 
2 is obtainable os a function of Vf, Z 0 , Z and latitude. 




(28) 




(29) 
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Fig. 4 

Geostrophic drag coefficient under neutral conditions as 
a function of surface Rossby number. 
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7. 



A revised formula for C/C a 

The revised formula for C a will be used in this section 
to compute C/C a , still usinp the previously tabulated values 
of peostrophic v/inds and C (table 3)« The results are shown 
in table 4 and tie resulting values cf C/C a then plotted 
versus Ri].4c as shown in f ip . 5. The procedure of fittinp 
the curve is shown below. 

After plotting the values of C/c a , a further restriction 
is placed on the empirical function 0/C a ; namely that C/C a sl 
at Ri=0. Moreover, since the shapes of the curves for C/C a 
in unstable and stable conditions are different laccordinp 
to the scatter of points), it is assumed that in stable condi- 
tions the desired function has the form y-s. e 9X where ys.C/C a 

and x ; iii^o* In unstable condition, y is assured to he of 

k x, 

the form y- l + k t X/£ * where y = C/C a and x = -Hi 140 » a ? k}, 
kQ, ar® constants to be deterr ined. 

(i) For stable conditions (table 5), with Rii40= x and 
C/C a =f y , if the assumed function has the form 



Bv the method of loast-s cuares , th normal eouation will be 



(30) 



then 



JL jy ^ a,* 



and if JLy* Y 



V - Ct)C 




( 31 ) 



Solving for a. 



a — Z|.i5^ 
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Substituting this constant in (30) and chan inf y to C/C a 
x to Hi]_ 4 o, the expression for C/C a Is obtained as 



C/Ca- ex.jp [- A 



(32) 
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TABLE ii. - Computations of C/C a using equation (28) for C a 



Bate CST 


log Ro 


r - 




Ca 


c 


6/Ca 


Ri il)0 


10 July 1956 


1305 


6.799 


-2.51 


0.081 


0.031)8 


0.0363 


1.01-5 


-0.1)6 


11)05 


6.773 


-2.50 


0.082 


0.0350 


0.0560 


1.600 


—0 0 17 


1505 


6.71)5 


-2.1)9 


0.083 


0.0360 


0.0521 


1. 1)1)5 


-0.28 


1530 


6.860 


-2.53 


0.080 


0.031)2 


0.0396 


1.160 


-0.25 


16o5 


6.730 


-2.1,8 


0.081) 


0.0359 


0.0530 


1.1)75 


-0.20 


1705 


6.699 


-2.1i7 


0.085 


0.0363 


0.061)8 


1.785 


-0.05 


1805 


6.687 


-2.1)6 


0.08.5 


0.0363 


0.0838 


2.320 


-0.02 


2105 


6.7U5 


-2.h9 


0.083 


0.0360 


0.0350 


0.975 


0.06 


23 July 1556 


0930 


6.500 


-2.1)0 


0.091 


0.0388 


0.0507 


1.305 


-0.15 


1105 


6.1)87 


-2.39 


0.092 


0.0393 


o.ob 5 i 


1.11)5 


-0.21 


1205 


6. 116I4 


-2.38 


0.093 


0.0397 


0.0515 


1.305 


-0.16 


1305 


6.1)39 


-2.37 


0.093 


0.0397 


0.0623 


1.570 


-0.09 


1605 


6.hl8 


-2.36 


0.091) 


0.01)01 


0.061)9 


1.615 


-o.li) 


1705 


6.hl3 


-2.36 


0.091) 


0.01)01 


0.071)6 


1.860 


-0.06 


1805 


6.107 


-2.36 


0.091) 


0.01)01 


0.0853 


2.150 


-0.02 


2130 


6.230 


-2.30 


0.100 


0.01)27 


0.0591 


1.3 pc 


O.Ol) 


2305 


6.0PI) 


-2.22 


0.109 


o.ol)65 


0.0850 


1.830 


0.03 


25 July 1956 


0205 


6.215 


-2.29 


0.101 


0.01)31 


0.0090 


0.209 


0.10 


oUo5 


6.516 


-2.1)0 


0.091 


0.0388 


0.0261) 


0.680 


o.o5 


0805 


6.721 


-2.1)7 


0.085 


0.0363 


0.01)21 


1.160 


-0.10 


0905 


6.9U8 


-2.56 


0.077 


0.0329 


0.0i)02 


1.220 


-0.01 


1230 


7.366 


-2.72 


0.066 


0.0282 


0.0368 


1.305 


-0.02 


il)05 


7.523 


- 2.11 


0.063 


0.0269 


0.0281) 


1.055 


-0.02 


1605 


7.593 


-2.80 


0.061 


0.0260 


0.0233 


0.895 


-0.01 


2330 


7.512 


- 2.11 


0.063 


0.0269 


0.0218 


0.810 


0.00 


6 Aug 1956 


1805 


7.188 


-2.65 


0.071 


0.0303 


0.0557 


1.835 


0.00 


2105 


7.196 


-2.66 


0.070 


0.0299 


0.0216 


0.722 


0.07 


2205 


7.201) 


-2.66 


0.070 


0.0299 


0.01)51) 


1.515 


0.01 


7 Aug 1956 


0U05 


6.295 


-2.32 


0.098 


0.0l)l8 


0.0970 


2.320 


0.01 


0805 


6.330 


-2.33 


0.097 


0.01)15 


0.0979 


2.360 


-o.ol) 


0905 


6.337 


-2.31) 


0.096 


0. 01)10 


0.1666 


i).06o 


-0.02 


1105 


6.350 


-2 .31) 


0.096 


0.01)10 


0.1032 


2.510 


-0.08 


1330 


7.111 


-2.62 


0.073 


0.0312 


0.0755 


2.1)20 


-o.ol) 


1605 


7.087 


-2.61 


0.071) 


0.0316 


0.08 61) 


2.730 


-0.01 


1805 


7.108 


-2.65 


0.071 


0.0306 


0.01)37 


1.51)0 


-0.01 


2005 


7.057 


-2.60 


0.071) 


0.0316 


0.031)1 


1.069 


o.o5 


2205 


7.000 


-2.56 


0.071) 


0.0329 


0.0153 


0.1)65 


0.13 
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TABLE lu (Cont'd) 



P= 



Date CST 


log fto 




e p 




C c/c a 


Ri lii0 


27 Aug 1956 


1230 


7.390 


-2.73 


0.065 


0.0278 


0.0325 1.170 


-0.07 


1505 


7.361 


-2.72 


0.066 


0.0282 


0.0307 1.090 


-0.02 


1705 


7.368 


-2.72 


0.066 


0.0282 


0.0318 1.125 


-0.02 


1930 


7.166 


—2 .61; 


0.071 


0.0303 


0.0023 0.076 


0.15 


2205 


7.392 


-2.72 


0.066 


0.0282 


0.0022 0.078 


0.0U 


29 Aug 1956 


2135 


6.377 


-2.35 


0.095 


O.OI1O6 


O.O669 1.6k) 


0.03 







TABLE 5. 


- C/C a in 


stable conditions 


• 

L 




c/Ca^y; 


S^y,- -V? 






1 


0.057 


0.975 


-0.025 


-0.001 


O.OC3 


2 


0.036 


1.385 


0.326 


0.012 


0.001 


3. 


0.035 


1.830 


0.601. 


0.021 


0.001 


U . 


0.099 


0.209 


-1.565 


-0.155 


0.010 


5. 


0.051 


0.680 


-0.386 


-0.020 


0.003 


6. 


0.067 


0.722 


-0.326 


-0.022 


o.ool. 


7. 


O.Olii 


1.515 


0.1.15 


0.006 


0.000 


8. 


0.008 


2.320 


0.81.2 


0.007 


0.000 


9. 


0.05I1 


1.069 


0.068 


o.ool. 


0.003 


10. 


0.127 


0.1-65 


-0.766 


-0.010 


0.016 


11. 


0.1*3 


0.076 


-0.27U 


-0.01.2 


0.023 


12. 


o.ol*o 


0.078 


-0.21*8 


-0.010 


0.002 


13. 


0.031 


1 .6I4O 


0.1-95 


0.015 


0.001 


urns 








-0.283 


0.068 
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(ii) For unstable conditions (table 6), the scatter 
of points of C/C a versus Rii40 suggests a curve of the 
form 

y=>+k,x,e kiX (3.3) 

where - Rii4C = x and C/C a = y. Again C/C a has been forced 
to approach unity as Ri approaches zero, If we let y-l-'f 
and then take natural logs 

,2m. ^ 2m. k, + 2 m x. +■ k.x- (34) 

= jlv, ♦ u 

)C 1 * 

Let 1L| = Y and 2vu k, s , and substitute in (34) 

Y - A. 4 k a X, 

The normal eouations are 





Y - a - - 0 


fW 

2Y*;- 

i-ci 


(W ^ 

0 

isi 1 * ' 



Solving for and k 2 

k 1= 30.364 , k 2 = -13. 324 

Substituting these constants in (33) and replacing x by 
-33-140, Y ty C /C a , gives 

C /Ce, - I - 30.36/, R lAo e)Cp [l3.32Z,TCi l/j0 ] (36) 

Note that C/C a of equation (36) may be maximized for Ri]_40 
- ' - -0.075, for which (C/C a ) max is 1.930. This is 

close to the maximum of the fitted curve of fig. 5. Lettau’s 
graph of C/C a has teen included in fig. 5 and it may be noted 

that the computations of t’'io paper are in general agreement 

22 



with those of Lettau, although his data did not span as 
a range of iii values. The results indicate that C/O o y 1 
stable conditions, whereas C/C a < 1 in stable conditions 



large 
in un- 
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l 

1 

2 

3 

a 

5 

6 

7 

f 

9 

10 

n 

12 

13 

i a 

15 

l6 

17 

18 

19 

20 

21 

22 

23 

2a 

25 

26 

Su-TT) 



o.a63 

0.165 

0.280 

0.251 

0.199 

0.053 

0.017 

0.ia7 

0.210 

0.]^7 

0.091 

o.iai 

o .058 
0.022 
0.103 
0.011 
0.017 
0.019 
o.oa2 
0.077 
o.oa3 
o.oil 
C.010 
0.075 

0.02a 

0.02a 

2.710 



TABLE 6 - C/C a in Unstable Conditions 



y.-=c/c^ 




ff/x,- 


i.oa5 


o.oa5 


0.097 


1.600 


0.600 


3.630 


i.aa5 


o.aa5 


1.590 


1.160 


0.160 


0.638 


i.a75 


o.a75 


2.390 


1.785 


0.785 


ia.700 


2.320 


1.320 


77.600 


1.305 


0.305 


2.080 


i.ia5 


o.ia5 


0.690 


1.305 


0.305 


i.9ao 


1,570 


0.570 


62.600 


1.615 


0.615 


a. 360 


1.860 


O.860 


ia.8oo 


2.150 


1.150 


52.310 


I.160 


0.160 


1.550 


1.220 


0.220 


20.000 


1.305 


0.305 


17.950 


1.055 


o.o55 


2.890 


2.360 


1.360 


32.aoo 


2.510 


1.510 


19.600 


2.a20 


i.a2o 


33.000 


2.730 


1.730 


157.300 


i.5ao 


o.5ao 


5a. 000 


1.170 


0.170 


22.700 


1.090 


0.090 


3.750 


1.125 


0.125 


5.200 







Yix; 


-2.333 


0.2ia 


-I.080 


1.289 


0.027 


0.213 


o.a 6 a 


0.078 


0.130 


-o.aa9 


0,063 


-0.113 


0.871 


o.oao 


0.173 


2.688 


0.003 


0.ia2 


a. 352 


0 


0.07a 


0.737 


0.022 


0.108 


-0.371 


o.oaa 


- 0.078 


0.663 


0.025 


0.10a 


a. 137 


0.008 


0.376 


1.172 


0.020 


0.208 


2.695 


0.003 


0.156 


3.957 


0 


0.087 


o.a38 


0.011 


o.oa 5 


2.996 


0 


0.033 


2.888 


0 


0.0 b 9 


1.061 


0 


0.020 


3.a78 


0.002 


0.1U6 


2.976 


0.006 


0.229 


3.500 


0.002 


0.151 


5.056 


0 


0.056 


3.989 


0 


o.oao 


3.122 


0.006 


0.23a 


1.322 


0.001 


0.032 


1.6a9 


0.001 


o.oao 


52.6a9 


0.576 


1.575 
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Pig. 5. 



Ratio c/c a as a function of Richardson number at Z = 140 cm, 
at O’Neill, Nebraska 1S56 with C a given by equation (28). 




8. The normalized — logarithmic wind sheai S 

It was explained in section 4 that if one knows 

S 3s — ^ - , one may obtain directly Cf/C s or 

O-* ' 

Ellison [2] has given a prediction equation for S in unstable 
cases 



5% ^ 3? = » 

where Z is height and L.' is the gradient length defined 



(37) 

as 



with 



y 



defined as 




q,e aG-/di 

kjj B© /dZ 



(58) 



r= r <39) 

here T is a constant, and and are eddy diffusivity 
coefficients for heat and momentum respectively 0 Dividing 
both sides of (37) by S 4 



Substituting 



(20) 



and 




= 5 



-4 



the expression for 



tz 



ana 



solving for S 



s = (■-■*'*.•) 



(40) 



Elackadar et al [2] suggest the value 18 for unstable 

conditions. However in stable conditions, they find that 
equation (40) £ ives no tetter approximation than the Lonin- 
Obukhov " log +• linear" wind profile defined by 



S - l + A 2 <X - constat* 

u ~ L 1 y (4i) 

which is ordinarily valid for Richardson numbers near cne c 
Again, it may be verified that (41) is equivalent to an 
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eouation 



5- l 4-c*’ (42) 
Blackadar et si have tabulated values of <*'' at C’lelll 

for various layers. For the layer 1-2 m under stable con- 
ditions, the best-fitting value of the constant was o<’~ ^ 
which will be used for the stable cases in this paper 0 

The computations for S are shown in table 7. The values 
of S are presumed applicable at the heirht assigned to Ri, 
that is at Z - Vi to or Z = 141.4 cm 0 The curve of S is 
shown in fig. 6. The result indicates that S (1 in unstable 
cases, and S > 1 in stable cases. Whether S in stable condi- 
tions levels off with a further increase of Ri, is not cer- 
tain yet, since the value chosen for =4* may not be valid for 
large positive values of Ri. 
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Results of computations for S 



TABLE 7. - 



Date CST 


^140 


Vg (kt s'] 


S ( 140) 


10 July 1956 


1305 


-0.46 


13.60 


0.572 


1405 


-0.17 


12.80 


0.708 


1505 


-0.28 


12.00 


0.638 


1530 


-0.25 


15.64 


0.652 


1605 


-0.20 


11.60 


0.683 


1705 


-0.05 


10.80 


0.850 


1805 


-0.02 


10.52 


0.935 


2105 


0.06 


12.00 


1.515 


23 July 1956 


0930 


-0.15 


19.44 


0.724 


1105 


-0.21 


19.00 


0.676 


1205 


-0.16 


18.00 


0.715 


1305 


-0.09 


17.00 


0 784 


1605 


-0.14 


16.20 


0.729 


1705 


-0.06 


16.00 


0.836 


1805 


-0.02 


15.80 


0.916 


2130 


0.04 


9.52 


1.326 


S305 


0.03 


7.52 


1.313 


25 July 1956 


0205 


0.10 


4.80 


1.888 


0405 


0.05 


9.60 


1.457 


0805 


-0.10 


24.00 


0.769 


0905 


-0.01 


26 . 00 


0.958 


1230 


-0.02 


32.00 


0.936 


1405 


-0.02 


46.00 


0.929 


1605 


-0.01 


54.00 


0.9^0 


2330 


0.00 


44.80 


1.033 


6 Aut' 1956 


1805 


0.00 


10 . 60 


0.994 


2105 


0.07 


10.80 


1.605 


2205 


0.01 


11.00 


1.123 


7 Aur 1956 


0405 


0.01 


11.88 


1.076 


0805 


-0.04 


12.88 


0.868 


0905 


-0.02 


13.08 


0.939 


1105 


-0.08 


13.48 


0.804 


1330 


-0.04 


15.84 


0.867 


1605 


-0.01 


15.00 


0.958 


1605 


-0.01 


15.76 


0.960 


2005 


0.05 


14.00 


1. 490 


2205 


0.13 


12.40 


2.147 
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TABLE 7. (Cont'd) 



Date CST 


^-140 


Vf (kts) 


S ( 140) 


27 Auf 1956 








1230 


-0.07 


26.48 


0.808 


1505 


-0.02 


24.80 


0.915 


1705 


-0.02 


25.20 


0.914 


1930 


0.15 


15.92 


2.375 


2205 


0.04 


26.60 


1.362 


29 Au f 1956 








2135 


0.03 


7.20 


1.278 
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The normalized logarithmic wind shear S as a function of 
Richardson number at Z=140 cm based on 1956 O'Neill, 
Nebraska data. 
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9. An alternative method of determination of U-* 

In section 7, C/C a v:as computed by eouations (8), (12), 
and (28). While there may be considerable error in assuming 
a functional form for C a of the type derived in section 6, 
there is also a possibility of error in the computation of 

by equation (8). Hence an alternative computation of 0*.* 
was employed, one which employs the expressions for S of the 
section 8. This approach essentially leads to a value of U. 4 
in terms of the wind-prof ile U.C5) and the inter ral-mean value 
of S in a sub-layer of the surface layer. Thus for example 

kz c>UL k 

~ U.*. 3 £ ~ () .W 2. 

and intepratinp from Z = Z 0 to Z - Z, assuming U-* is constant 
in the surface layer leads to 




where 3 is the intep ral-mean value of S up to level Z. 

As in the previous section S was taken in accordance 

with 

3 \ — l8 unstable cases 

( 44 ) 

■5— * "** 9 stable cases 

Although Priestly |[8] has raised a euestion re; ardinp the 
the validity of a constant U.* with hei ht in stable condi=- 
tions, it was assumed that a surface layer at the depth of 
at least 2m exists for both stable and unstable cases. 

Values of S accordinp to eouation (44) were plotted 
usinp Hi (Z) for values of Z- 5, 35, 70, 140, 280cns. The 
lowest value of heipht corresponds to a level at the eonetric 
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mean of the roughness parameter Z 0 and Z = 25 cm, usir.p the 
mown values of and do/dz for this layer,, The 

praphs of S versus ^ Z was plotted and the integral-mean 
value S extracted. Fip. 7 shows an example of an unstable 
and stable case which occurred at 1905 and 2105 CST, respec- 
tively on 10 July 1956. This permits one to solve for b\ 
means of 

g U(iQ.) U5) 

S Jwv (a©o/z 0 ) 

The values of U-* computed by eouation (45) are displayed in 
table 8 alonf with the values of U**, listed earlier in 
table 1. The absolute averare error is approximately 10$, 
although the alpebraic error is 3$. 

Note that the ratio Cf/c a may be obtained as the product 

C£ _ C£ c_ 

c«. ■ C C 4 

— Cj ^5 (45) 

The value of S as a function of fti(Z) may be obtained from 
eouation (44). Values of m* may be obtained from equation 
(45); and C a may be determined from eouation (28) „ 

The main purpose of this last discussion is to empha- 
size that Cf is not eoual to C a . Listed below in table 8 
are some computed valu.es of Cf/C a „ 
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The normalized logarithmic wind shear S as a function of 

for unstable and stable conditions for 10 July 1S56. 



TABLE 8 



C f /C a in unstable and stable cases 



Date CST 


l?i 140 


S 


u £: i / sec ) 


(Cf/’C a ) U^rr/sec) f rom ( 


10 July 1956 


1305 


-0.46 


0.572 


25.1 


0.605 


24.7 


1405 


-0.17 


0.9C0 


35.0 


1.105 


Q 

ksk) o *• 


.505 


-0.28 


0.885 


31.3 


0.925 


31.2 


1605 


- 0.20 


0 . 900 


31.6 


1.038 


35.8 


1705 


-0.05 


0.950 


34.8 


1.515 


35.0 


1805 


-0.02 


0.97C 


43.0 


2.100 


44.0 


2] 05 


0.06 


1.150 


20.5 


1.435 


21.0 


23 July 1956 


1105 


- 0.21 


0.850 


35.7 


0.648 


42.9 


1205 


-0.16 


0.880 


38.7 


0.773 


46.7 


1305 


-0.09 


0.685 


44.3 


1.025 


52.9 


1605 


-0.24 


0.885 


42.6 


0.956 


52.6 


1705 


-0.06 


0.950 


46.7 


1.215 


59.7 


2305 


0.04 


1.115 


27.1 


2.030 


31.8 


25 July 1956 


0405 


0.05 


1.100 


16.9 


1.325 


12.7 


6 Aup 1956 


2105 


0.07 


1.200 


12.1 


1.210 


11.7 


7 Aup 1955 


2005 


0.05 


1. 100 


22.6 


1.517 


27 o 


27 Aur 1956 


2205 


0.04 


1.00 


10.6 


0.386 


7.98 
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10. Summary and conclusions. 

A functional relationship betv/een C/C a and Ri was 
obtained by using an expression for C a suggested by Lettau. 
Better agreement was obtained after making some revisions 
in Lettau's formula for C a . An expression for C a has the 
practical significance that one can pet directly the ratio 
of a windspeed at some anemometer height to a surface geo- 
strophic wind by knowing the roughness parameter. 

By using the revised expression for C a , a new functional 
relationship between C/C a and Ri was obtained. This function 
has different analytical forms in unstable and stable condi- 
tions. 

Another way of obtaining C/C a was based upon computa- 
tion of U.* using Blackadar’s results. By this means a 
wind-profile expression was obtained which varies with 
stability. Prom this new procedure Dne can ret the value 
of Cf/C a without using C/C a . 

s 
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APPENDIX 



Computation for C a usinp Lettau’s formula 



The original ecuation is 



C,, - 



o.'oA 



( Ca, R«») — 

By cross multiplying 

Ioq(C^o) ^ + o.voa/c«. 



or 



Cc^Ro = »o 



(2., 24 +0 »c4/Ca) 



OUC4/C t 



th en 



Co.R 



73>S 



s 173 ,^ X iO 



/ 0 . 10^1 n i.3o3 \ 

° = ex P ( ; 



“ exp ( O-i’^’j/Ca) 

Let 0 o 2395/Ca^y and substi tutinp in the above eolation 

fcv) = -Bs- - e y ~ o 

" 73.8 \y/ 

__ Q, 23 fp>/ 173. g 'Ro - V g* _ q 

3/ 



Since y ^ O 



fly) - o»i3ss/i]3;8R 0 -ye*-o (4?) 

Takinr the lops to the base 10 

^ 3 (££i22'R„) = L.,y + 0-13437 (48) 

Let lop l< — constant, and replace lop y by 

0 e 4242 jU-y and substitute in (48) 

Ay) — o ^ y + o-4^iy - i< = © 

To pet the root of ecuation (49), we will use the 



( 49 ) 



Newton- riaphson method 



y\u = v f - 



f(y-) 



(50) 



Since f Cy)- 0.4343/y4- 0.4343, therefore by substituting this 
in (50) we obtain 

y,. - y. _ Q>M343jUy? 4- - K 

Substituting lop Roj for K and divide the fraction 

by C. 4343 , 

0.4343 

Vu, = y,- - Av.- » y. - I Am h 3 L„ 

Chanpinp 2wy. to 2,3026 lop yj_ 



y = y - 

•'l + i y i 



6 y," + y,- 4 - 6* $3^.7 ~ R© 



»/y,* +• » 



(51) 
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